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The Generation of Magnetic Fields and X-ray Observations 
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Abstract. We show that strong magnetic fields can be generated at shock waves associated with formation of galaxies or 
clusters of galaxies by the Weibel instability, an instability in collisionless plasmas. The estimated strength of the magnetic 
field generated through this mechanism is close to the order of values observed in galaxies or clusters of galaxies at present, 
which indicates that strong amplification of magnetic fields after formation of galaxies or clusters of galaxies is not required. 
This mechanism could have worked even at a redshift of ~ 10, and therefore the generated magnetic fields may have affected 
the formation of stars at the early universe. This model will be confirmed by future observations of nearby clusters of galaxies. 
In this context, we also present the Japanese X-ray missions. 
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1. Introduction 

The question of the origin of galactic magnetic fields is one of 
the most challenging problems in modern astrophysics. It is 
generally assumed that magnetic fields in spiral galaxies are 
amplified and maintained by a dynamo through rotation of 
the galaxies (Widrow 2002). The dynamo requires seed fields 
to be amplified. However, observations of microgauss fields 
in galaxies at moderate redshifts strongly constrain the lower 
boundary of the seed fields (Athreya et al. 1998). Moreover, 
magnetic fields are also observed in elliptical galaxies and 
galaxy clusters where rotation, and thus dynamo mechanism 
cannot play a central role (Clarke, Kronberg & Bohringer 
2001; Widrow 2002; Vallee 2004). 

On the other hand, the Weibel instability can provide an- 
other mechanism to generate strong magnetic fields (Weibel 
1959; Fried 1959). This instability is driven in a collisionless 
plasma, or a tenuous ionized gas, by the anisotropy of the 
particle velocity distribution function (PDF) of the plasma. 
When the PDF is anisotropic, currents and then magnetic 
fields are generated in the plasma so that the plasma parti- 
cles are deflected and the PDF becomes isotropic (Medvedev 
& Loeb 1999). Through the instability, the free energy at- 
tributed to the PDF anisotropy is transferred to magnetic field 
energy. This instability does not need seed magnetic fields. 
It can be saturated only by nonlinear effects, and thus the 



magnetic fields can be amplified to very high values. This 
instability has been observed directly in recent laser experi- 
ments (Wei et al. 2002). In astrophysical plasmas, the insta- 
bility is expected to develop at shocks or at steep temper- 
ature gradients, where the PDF is anisotropic. Examples of 
the sites are pulsar winds, shocks produced by gamma-ray 
bursts, jets from active galactic nuclei (AGNs), cosmological 
shocks, and cold fronts (contact discontinuities between cold 
and hot gas) in clusters of galaxies (Medvedev & Loeb 1999; 
Kazimura et al. 1998; Nishikawa et al. 2003; Schlickeiser & 
Shukla 2003; Okabe & Hattori 2003). Although the instabil- 
ity was found in 1959, its nonlinear nature had prevented us 
from understanding its long-term evolution. Recently, how- 
ever, as computer power increases, detailed particle simula- 
tions of plasmas have been initiated and they have revealed 
the evolution of magnetic fields even after saturation of the 
instability (Silva et al. 2003; Medvedev et al. 2005). Based on 
these results, we consider the generation of magnetic fields at 
galaxy and cluster-scale shocks through the Weibel instability 
at the formation of galaxies (both ellipticals and spirals) and 
clusters. We use the cosmological parameters of Hq = 0.3, 
A = 0.7, the Hubble constant of Hq = 70 km s"^ Mpc"\ 
and (jg = 0.9. 
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2 MODELS 



2.2 Shock formation 



2. Models 

2.1. Generation of magnetic fields 

The evolution of the Weibel instability can be understood 
as evolution of current filaments generated in a plasma 
(Medvedev et al. 2005; Kato 2005). The growth of the insta- 
biUty ceases once the magnitude of current in each filament 
reaches the Alfven current (Kato 2005); physically this con- 
dition is consistent with the condition by Medvedev & Loeb 
(1999) that the particle's gyroradii in the excited magnetic 
fields are comparable to the characteristic wavelength of the 
magnetic field. In proton-electron plasmas, heavier protons 
mainly determine the final magnetic field strength (Frederik- 
sen et al. 2004). Thus, we consider the generation of magnetic 
fields by protons. At a shock front, an anisotropy is produced 
in the PDF of protons owing to the difference between the 
shock speed u and the thermal velocity of protons vth- When 
the shock Mach number (M = u/vth) is larger than ~ 2, 
which is the case in the galactic shocks we consider later, the 
anisotropy becomes large enough to adopt the 'particle-limit' 
solution discussed by Kato (2005). In this case, the saturated 
magnetic field is determined only by the shock velocity and 
the proton number density Hp as 

-Bsat ~ XPU^f'^T^npTTlp , (1) 

where rup is the proton mass, and Xp is the isotropization fac- 
tor at the saturation, which is typically 0.5 (Kato 2005). The 
saturated magnetic energy density reaches a sub-equipartition 
level with the particle kinetic energy, which is consistent with 
previous simulation results under large anisotropic conditions 
(Yang, Arons & Langdon 1994; Califano et al. 1998) 

After saturation, the magnetic field strength decreases 
as the current filaments merge. Recent particle simulations 
have shown that the final magnetic field strength is given 

1 /2 

by Bf = r/mor-Bsat, where rj^ci ^ 0.01 (Silva et al. 2003; 
Medvedev et al. 2005), thus we use this value. Long-term 
evolution (say Gyr) of the generated magnetic fields is stiU 
an open question. However, Medvedev et al. (2005) indicated 
that the fields do not decay rapidly after the saturation, be- 
cause the spatial scale of the fields increases rapidly and dis- 
sipative processes such as Ohmic heating or ion-neutral col- 
lisions, if any, are suppressed. The very long-term evolution 
of the magnetic fields, for which the results of Medvedev et 
al. (2005) and Silva et al. (2003) cannot be directiy applied, 
is discussed in Fujita and Kato (2005). 

2.2. Sliock formation 

According to the standard hierarchical clustering scenario of 

the universe, an initial density fluctuation of dark matter in 
the universe gravitationally grows and collapses; its evolu- 
tion can be approximated by that of a spherical uniform over- 
dense region (Gunn & Gott 1972; Peebles 1980). The col- 
lapsed objects are called 'dark halos' and the gas in these 
objects later forms galaxies or clusters of galaxies. At the 
collapse, the gas is heated by the 'virial shocks' to the virial 
temperature of the dark halo, Tvir = GM/ (2rvir), where G 
is the gravitational constant, and M and Tvir are the mass and 



the virial radius of the dark halo, respectively. The relation 
between the virial radius and the virial mass of an object is 
given by 



^vir — 



3M 



1/3 



(2) 



where Pcrit(2^) is the critical density of the universe and 
Ac(z) is the ratio of the average density of the object to the 
critical density at redshift z. In the ACDM cosmology, the 
critical density depends on redshift because the Hubble con- 
stant depends on that, and it is given by 



Pcrit,0^o(l + Z)^ 

n{z) 



(3) 



where pcrit.o is the critical density at z = 0, and Q{z) is the 
cosmological density parameter given by 



n{z) 



no{l + zf 



(4) 



for the flat universe with non-zero cosmological constant. 
The ratio Ac{z) is given by 



Ac{z) = 18 TT^ -I- 82a; - 39a;^ 



(5) 



for the flat universe (Bryan & Norman 1998), where the pa- 
rameter X is given by a; = — 1. The virial shocks form 
at r w Tvir and the velocity is Wvir ~ -\/GM/Vvir. 

In addition, recent cosmological numerical simulations 
have shown that 'large-scale structure (LSS) shocks' form 
even before the gravitational collapse (Miniati et al. 2000; 
Ryu et al. 2003; Furlanetto & Loeb 2004); they form at flie 
turnaround radius (rta ~ 2rvir), the point at which the den- 
sity fluctuation breaks off from the cosmological expansion. 
Thus, the gas that is going to form a galaxy or a cluster under- 
goes two types of shocks; first, the gas passes the outer LSS 
shock, and then, the inner virial shock. The typical velocity of 
the LSS shocks is w H{z)Rp (Furlanetto & Loeb 2004), 
where H{z) is the Hubble constant at redshift z, and Rp is 
the physical radius that the region would have had if it had 
expanded uniformly with the cosmological expansion. The 
temperature of the postshock gas is Tg « 3/16 {nmp/kB)vg, 
where fim,p is the mean particle mass, and ks is the Boltz- 
mann constant. Note that we do not consider mergers of ob- 
jects that have already collapsed as the sites of magnetic field 
generation. During the merger, collapsed objects just bring 
their magnetic fields to the newly bom merged object. 

Since the Weibel instability develops in ionized gas 
(plasma), we need to consider the ionization history of the 
universe. After the entire universe is ionized by stars and/or 
AGNs (z ^ 8), magnetic fields are first generated at the LSS 
shocks. In this case, we do not consider the subsequent gener- 
ation of magnetic fields at the inner virial shocks, because the 
strength is at most comparable to that of the magnetic fields 
generated at the LSS shocks. On the other hand, when the 
universe is not ionized (z ^ 8), the Weibel instability cannot 
develop at the outer LSS shocks. However, if the LSS shocks 
heat the gas (mostly hydrogen) to > 10^ K and ionize it, 
the instability can develop at the iimer virial shocks. 



3 RES ULTS AND DISCUSSION 



For this case, the gas ionized at the LSS shock may re- 
combine before it reaches the virial shock. The recombina- 
tion time-scale is given by 



0.7 / N 



,(6) 
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r,ec = «1.22xl0^yr- — - 
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where a is the recombination coefficient, is the electron 
density, T is the gas temperature, y is the ionization frac- 
tion, and uh is the hydrogen density (Shapiro & Kang 1987). 
If we assume that Ticc — Tdyn, where Tdyn ~ (l/2)rta/u 
is the time-scale that the gas moves from the LSS shock to 
the virial shock, the ionization rate when the gas reaches the 
virial shock is 



y 



^dyn 



1.22 X 105 yr 



VIO^K J 



(7) 



for y < I. We found that y < 1 for z ^ 9, and the mini- 
mum value when the generation of magnetic fields is effec- 
tive (z ^ 12, see ^ is y ^ 0.3. When j/ < 1, we simply 
replace Up in equation Q with yuH- For temperature, we as- 
sumed T = Ts in equation Q. Since the magnetic fields do 
not much depend on the temperature (Bgat oc y'~'^^ oc T^^^), 
they do not much change even when radiative cooling reduces 
the temperature; at least i?sat does not change by many or- 
ders of magnitude. Since the ionization rate is fairly large, the 
ambipolar diffusion of magnetic fields can be ignored (e.g. 
eq. 13-57 in Spitzer 1978). It was shown that the ionization 
rate just behind a shock is y ~ 0.1, if the shock velocity is 
relatively small 40 km s^^ (Shapiro & Kang 1987; Susa et 
al. 1998) . If the shock velocity is larger, j/ ^ 1. In our cal- 
culations, the velocity of the LSS shocks is u ^ 40 km s~^, 
when the generation of magnetic fields is effective. Thus, y 
just behind the shocks is at least comparable to that obtained 
through the condition of Tree ^ Tdyn, and the incomplete ion- 
ization does not much affect the results shown in the next 
section. 



3. Results and discussion 

FigQ shows the typical mass of objects, M, as a function 
of redshift z; the labels la, 2a, and 3cr indicate the am- 
plitudes of initial density fluctuations in the universe from 
which the objects form, on the assumption of the CDM fluc- 
tuations spectrum (Barkana & Loeb 2001). Fig|2]shows the 
downstream temperature at the virial shock (Tvii) and that 
at the LSS shock (Ts) for the objects; Tvir is always larger 
than Ts- The ratio T^n/Ts indicates that > 4 for the virial 
shock. For the LSS shock, ^ 1, because the gas out- 
side the shock is cold. Thus, equation Q can be applied to 
both shocks. In Fig|3l we present the strength of magnetic 
fields {Be) at a scale of r^-„ for the collapsed objects. We as- 
sume that the entire universe is reionized at z = 8. Thus, for 
z > 8, magnetic fields are generated only at the virial shocks 
if Ts > 10"* K. We assume that Be ^ Bf and plot the lines 
only when Tg > 10* K. The recombination is effective at 
z > 12 for the 3a model and at z > 8 for the 2a model. On 
the other hand, for z < 8, the magnetic fields are generated 
at the LSS shocks. We consider the compression of the fields 
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Fig. 1. Typical masses of objects forming at redshift z. The 
labels, la, 2a, and 3a (solid, dotted, and dashed lines, re- 
spectively) indicate the amplitudes of initial density fluctu- 
ations from which the objects formed; ^ l-3a is the typ- 
ical value (Barkana & Loeb 2001). Objects with masses of 
M < 10^2 Mq and M > 10^^ M© correspond to galaxies 
(ellipticals and spirals) and clusters of galaxies, respectively. 



while the size of the gas sphere decreases from r = rta to 
r = Tvir, and thus we assume that Be = 8^^^ Bf. Moreover, 
we plot the lines only for Tvir > 2 x 10^ K, because below 
this temperature, gas infall is suppressed by photoionization 
heating (Efstathiou 1992; Furlanetto & Loeb 2004). 

In Figls] the strength of magnetic fields reaches ^ 10^^- 
lO"'^ G, which is very close to the values observed in nearby 
galaxies and clusters of galaxies (^ 10~^ G; Clarke et al. 
2001; Widrow 2002; Vallee 2004). This indicates that strong 
amplification of magnetic fields, such as dynamo amplifica- 
tion, is not required after formation of the galaxies and clus- 
ters. The absence of a strong amplification is consistent with 
the observations of galactic magnetic fields nt z^ 2 (Athreya 
et al. 1998). Future observations of higher-redshift galaxies 
would discriminate between our model and strong dynamo 
amplification models; the latter predict much weaker mag- 
netic fields at higher redshifts. Moreover, since the predicted 
galactic magnetic fields are comparable to those at present, 
they might have affected the formation of stars in protogalax- 
ies. Fig|3 also shows that our model naturally explains the 
observational fact that the magnetic field strengths of galax- 
ies and galaxy clusters fall in a small range (a factor of 10). 
Since our model predicts that magnetic fields are generated 
around objects, magnetic fields in intergalactic space are not 
required as the seed or origin of galactic magnetic fields. 

Although it would be difficult to directly observe the gen- 
eration of magnetic fields through the Weibel instability for 
distant high-redshift galaxies, it would be easier for nearby 
clusters of galaxies. Since clusters are now in the forma- 



4 JAPANESE X-RAY MISSIONS 



o 




1+z 

Fig. 2. Temperatures behind the virial shocks (Tviri thick 
hnes) and those behind the LSS shocks (Ts ; thin Hnes) for ob- 
jects forming atredshift z (Fujita & Kato 2005). SoHd, dotted, 
and dashed lines correspond to la, 2a, and Su fluctuations, 
respectively (see FigQJ 
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Fig. 3. Typical magnetic field strengths of objects forming at 
redshift z (Fujita & Kato 2005). Solid, dotted, and dashed 
lines correspond to \a, 2a, and 3a fluctuations, respectively 
(see Fig[l). 



tion process, LSS shocks should be developing outside of 
the virial radii of the clusters (Miniati et al. 2000; Ryu et 
al. 2003). Since particles are often accelerated at shocks, the 
synchrotron emission from the accelerated particles could 
be observed with radio telescopes with high sensitivity at 
low frequencies (Keshet, Waxman & Loeb 2004). The total 
non-thermal luminosity due to synchrotron luminosity plus 
inverse Compton scattering with cosmic microwave back- 
ground (CMB) photons is estimated as 



- in 



/Ml 2 

7:U 



rta/u mp 2 



(8) 



where e is the acceleration efficiency and / is the gas frac- 
tion of a cluster. If we assume e = 0.03 and / — 0.15, 
the maximum value of L^t is ~ 10*'^ erg s~^. Since the 
energy density of magnetic fields (ub) is smaller than that 
of the CMB (ucmb), most of the non-thermal luminosity 
(Lnt) is attributed to the inverse Compton scattering such as 
ucMBint/(uB + mcmb), which may have been detected in 
the hard X-ray band (^20 keV; Fusco-Femiano et al. 2004). 
Thus, the synchrotron radio luminosity is UBLnt/ {ub + 
ucmb) ^ 10^^ erg s^^. Some of the diffuse radio sources 
observed in the peripheral cluster regions ('radio relics') may 
be this emission (Govoni & Feretti 2004). Since the Weibel 
instability generates magnetic fields on the plane of the shock 
front, the synchrotron emission should be polarized perpen- 
dicular to the shock front (Medvedev & Loeb 1999), which 
is actually observed for some radio relics (Govoni & Feretti 
2004). The synchrotron emission will tell us the positions of 
the LSS shocks, if they exist. If magnetic fields are gener- 
ated there, they should be observed only downstream of the 
shock. This may be confirmed through Faraday rotation mea- 
surements of radio sources behind the cluster for both sides 
of the shock, if the coherent length of the fields sufficiently 
increases. 



4. Japanese X-ray missions 
4.1. Suzaku (Astro-E2) 

Suzaku is the recovery mission for ASTRO-E, which did not 
achieve orbit during the launch in 2000. Suzaku was launched 
in July 10, 2005. Suzaku has three detectors (XRS, XIS, and 
HXD). XRS has a superb energy resolution ( ^ 6.5 eV 
^ 100 km s^^ ). We expected that it could detect gas mo- 
tion in clusters, for example. XRS uses liquid helium to cool 
the detector. Unfortunately, the helium evaporated before the 
first light (Aug. 8) and XRS was lost. The cause of the trouble 
is now under investigation. The observation schedule will be 
changed and optimized to the remaining two detectors (XIS 
and HXD). XIS has imaging capability and its energy range 
is 0.2-12 keV. HXD is a hard X-ray detector and its energy 
range is 10-600 keV. XIS and HXD are working well. Their 
sensitives are shown in Figure 0] It is to be noted that the 
background levels of these detectors are lower than those of 
Chandra and XMM-Newton. XIS and HXD will give us vari- 
ous information in the field of high energy astrophysics (e.g. 
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hard X-ray from black holes and clusters of galaxies). Read- 
ers can get the information on Suzaku at several web sites.' 



Fig. 5. A simulated hard X-ray image and the spectrum of 
the radio relic in Abell 2256 observed by NeXT. 



4.2. NeXT 

NeXT (New X-ray Telescope ) is the major Japanese X-ray 
mission after Suzaku} If the plan is approved, NeXT will be 
launched in 2011, before the launch of Constellation-X and 
XEUS. NeXT will have imaging capability in the hard X-ray 
band. For example, we will easily get the images of radio 
relics in clusters of galaxies by observing their inverse Comp- 
ton emission (Figure|5jl. The sensitivity will be much higher 
than those of previous satellites (Figure |6j. NeXT will also 
have an ultra-high resolution soft X-ray spectrometer, which 
will be another challenge of Suzaku XRS. Moreover, NeXT 
will have a soft gamma-ray detector Thus, it will enable us 
to observe objects in a very wide energy band. 
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